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We prescribe a pseudo-Newtonian vector potential for studying accretion disks 
. around Kerr black holes. The potential is useful to study the inner properties 

of disk not confined to the equatorial plane where general relativistic effect is 
indispensable. Therefore, we incorporate the essential properties of the metric at 
the inner radii through the pseudo-Newtonian potential derived from the general 
Kerr spacetime. The potential, reproducing most of the salient features of the 
general-relativity, is valid for entire regime of Kerr parameter. It reproduces the 
last stable circular orbit exactly as that in the Kerr geometry. It also reproduces 
last bound orbit and energy at last stable circular orbit with a maximum error 
~ 7% and ~ 15% respectively upto an orbital inclination 30°. 
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1. Introduction 

The standard model of the accretion disk proposed by Shakura & Sunyaev (1973) ex- 
plains the dynamics of a thin disk by the Newtonian potential, not taking into account of 
the essential general-relativistic effects of the black hole. The only notion it hires from gen- 
eral relativity (hereinafter GR) is that the inner most edge of the disk truncates at the last 
stable circular orbit of the Schwarzschild geometry. Complete general relativistic descrip- 
tion of the thin Keplerian disk was developed by Novikov & Thorne (1973). Owing to the 
complicated and cumbersome nature of underlying general relativistic equations, especially 
to study the full scale dynamics of the disk, it would be worthy enough to study in Newto- 
nian framework using a modified Newtonian potential, namely pseudo-Newtonian potential 
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(hereinafter PNP), having all the essence of a black hole geometry. Such a potential was first 
proposed by Paczyhski & Wiita (1980; hereinafter PW) for modeling a thick/thin accretion 
disk around a nonrotating black hole. The potential exactly reproduces marginally stable 
orbit (r s ) and marginally bound orbit (r^) of that in full GR and the efficiency per unit 
mass at the last stable circular orbit (E s ) agrees with at most 10% error. Nowak & Wagoner 
(1991) and Mukhopadhyay & Misra (2003) proposed other potentials for an accretion disk 
to describe epicyclic frequency which mimic most of the properties of the disk governed by 
the Schwarzschild and Kerr geometry respectively. 

PNPs to describe the fluid dynamics of accretion disk around a rotating black hole 
in equatorial plane were also proposed by Chakrabarti & Khanna (1992) and Artemova et 
al. (1996) which, however, were not valid for entire range of Kerr parameter. Later, a 
twodimensional potential was described by Mukhopadhyay (2002; hereinafter M02) which is 
valid for the entire range and reproduces r s exactly, r& with a maximum error ~ 5% and E s 
with a maximum possible error ~ 10% compared to that in the Kerr geometry. The more 
interesting fact, however, lies in the methodology adopted by M02 in deriving the potential 
which can be used to formulate the PNP for any metric according to the physics concern. 
Based on the procedure, several authors have developed PNPs at the equatorial plane for 
other metric, e.g Hartle-Thorne metric by Ghosh (2004), Kerr- Newman metric by Ivanov 
& Prodanov (2005). We plan to follow the procedure outlined by M02 to derive a PNP 
for general Kerr geometry (particle motion not confined to the equatorial plane) which can 
be applied to study the generalized axisymmetric accretion disk dynamics. Naturally, the 
aim is to furnish a pseudo-Newtonian vector potential rather than a scalar potential with 
radial and vertical component in cylindrical coordinate system suitable for studying a thick 
accretion disk. 

While some authors focus on to propose an ad-hoc and a mathematically simpler form 
of PNP, we derive it from the spacetime metric directly. Therefore, we do not need to worry 
about its final form strictly determined by our procedure for any spacetime metric unlike the 
other PNPs. Problem with an ad-hoc PNP is that there is no easy one to one correspondence 
between spacetime metric and PNP. 

Most of the cosmic objects in our nature are rotating (see e.g. Iwasawa et al. 1996; 
Wilms et al. 2001; Miniutti et al. 2004; Miller et al. 2004; Genzel et al. 2003). Re- 
cently, it has also been argued (McClintock et al. 2006) that the source GRS 1915+105 is 
a rapidly rotating black hole with spin parameter a > 0.98. Thus, the predictions of the 
inner disk properties will be incorrect if the pseudo-Newtonian modeling fails to include the 
spin parameter of the black hole. The potential proposed by M02 has already been used 
to study the hydrodynamics of the disk around Kerr black holes, which shows a dramatic 
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change of the valid parameter space region; any shock location and other fluid properties 
change for different rotation parameters of the black hole (Mukhopadhyay 2003). Recently 
Chan et al. (2005, 2006) have used this potential for time dependent studies of accretion 
flow for two dimensional viscous hydrodynamic disk and for threedimensional self-gravitating 
magnet ohydro dynamic disk around rotating black holes using spectral method. 

Our aim in the present paper is to develop a potential for the general Kerr geometry 
which can be used to mimic the general relativistic effects of the spacetime, enabling us to 
suitably use it to study the accretion dynamics, especially of a thick disk with the inclusion 
of vertical outflow. Thus we goal for a generalized potential which is a function of r and 9 
both and having two components, in radial and vertical direction, in cylindrical coordinate 
system, hence a vector potential. This is useful for both quasi-numerical analyzes and time 
dependent simulations, avoiding the use of complicated relativistic equations, yet capturing 
the essential physics of spacetime outside of a Kerr black hole. In early, some attempt 
was made to describe PNP in generalized Kerr geometry by, e.g., Dadhich (1985), Semerak 
& Karas (1999), Chakrabarti & Mondal (2006). The PNPs described in first two works, 
while aim to furnish frame-dragging effects in the Newtonian framework, do not reproduce 
marginally stable and bound circular orbits and specific energy close to the black hole and 
hence may not be suitable for the hydrodynamical analysis of a Keplerian accretion disk, 
particularly, when aim is to study inner disk properties. The third one is the modification 
of Chakrabarti & Khanna (1992) potential to mimic approximately the general relativistic 
effects in accretion disks for — 1 < a < 0.8. The potential is a prescription with which the 
authors maintain the ad-hoc-ness of PW potential and may be used to study accretion disk 
having particle orbits of small inclination angle at a certain parameter window. Our goal 
here, however, is to extend the procedure initiated by M02 which can be used to describe a 
PNP for the general Kerr spacetime; particle motion not confined to the equatorial plane. 
The potential given by M02 is more useful to study the dynamics of a vertically averaged thin 
disk at equatorial plane. Henceforth, we wish to derive a generalized potential applicable 
for off- equatorial planes following the procedure of M02 which reproduces exactly or in good 
agreement all the (inner) accretion disk properties (unlike the previous works) for a general 
case in Kerr geometry. This will be useful to study the thick accretion disk dynamics with 
the inclusion of vertical outflow and wind. The potential should reproduce those features of a 
rotating black hole geometry which have been reproduced by that of M02 at equatorial plane 
and that of PW for a nonrotating black hole. We formulate our PNP from the Kerr metric. 
As the metric is involved directly to our calculation, many of the features of Kerr geometry 
are inherent in our potential. The axisymmetric PNP we plan to establish here for a general 
Kerr metric is more useful for quasi-analytical studies and time dependent simulations of 
accretion and wind and related physics as that studied by Chan et al. (2005, 2006), Lipunov 
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& Gorbovskoy (2007), Shafee et al. (2007) using M02, without loss of any generality. 

We arrange the paper in the following manner. In the next section, we present the 
basic equations and derive the PNP. In §3, we compare the various features of dynamics of 
accreting particles obtained by our potential with that in Kerr geometry. Finally §4 presents 
a discussion with summary. 



2. Formulation of the generalized pseudo-Newtonian potential 

Here we furnish the derivation of our said potential. It is condensed in three subsections, 
where we have elucidated the description of the potential starting from the Kerr metric. 



2.1. Integrals of motion of a test particle in the Kerr geometry 

The Kerr spacetime in the Boyer-Lindquist coordinate system in geometrized units 
G = c = 1 is given by 



/ 2Mr\ , 2 4aMrsin 2 # , ,, S l2 
ds = — f 1 — J dt 2 dtd(f) + —dr 2 



^ ,„9 ( 9 9 2Mra 2 sin 2 6\ 9 n , , 9 
+Zd6 2 + ir 2 + a 2 + Jsin 2 #d0 2 , (1) 

where A = r 2 + a 2 — 2Mr, E = r 2 + a 2 cos 2 9 and a = which is specific angular momentum 
of the black hole named Kerr parameter. The Lagrangian density of the particle of mass m 
in the Kerr spacetime is given by 

/ 2Mr\ . 2 4aMrsin 2 #. • S 2 
2C = -(1- —) t 2 10+ -f 2 



^^9/9 9 2Mra 2 sin 2 6»\ .2^2 
+W 2 +ir 2 + a 2 + Jsin 2 ^ 2 , (2) 

where over-dots denote the derivative with respect to proper-time r along the curve. From 
the symmetries, we obtain two constants of motion corresponding to two ignorable coordi- 
nates t and given by 

dC / 2Mr\ . 2aMrsin 2 # ; 
Pt = = - II — J t (f) = constant = —E (3) 

and 

dC 2aMrsin 2 6. ( 2 2 2Mra 2 sin 2 ^\ 2 n ■ 

= — - = — 1 + I r + a H — sin Ocp = constant = A 2 , (4) 

d(j) £ \ L / 



where X z is the component of angular momentum of the orbiting particle along spin axis of 
the black hole (i.e. parallel to the symmetry axis). The Hamiltonian of the system can also 
be written as 

n = \g» u p»Pv. (5) 

The Hamiltonian is itself a constant of motion and by usual setting g^vjf'p 1 ' = —m 2 for a 
particle with nonzero rest mass, we obtain the third constant of motion (integral of motion) 

as 

TC = -m 2 = £ = constant. (6) 

The Lagrangian and the Hamiltonian of the geodesic with the above condition are exactly 
same. The three integrals of motion obtained above are sufficient to determine the geodesic 
equations under restrictive condition like the motion of a test particle confined to an equa- 
torial plane {9 = 0, 9 = ir/2) (e.g. Shapiro & Teukolsky 1983). To describe the motion of 
a particle for a generic case around the Kerr black hole, with dependency on 9, a fourth 
integral of motion is needed which does not come from the symmetries of the metric. It 
was discovered by Carter (1968) that such an integral of motion can be obtained explicitly 
by separability of the Hamilton- Jacobi equation (hereinafter HJE) (also see Chandrasekhar 
1983) given by 

Q =p 2 e + cos 2 9 [a 2 (m 2 - E 2 ) + -^-}. (7) 

sin 9 

Q is called the Carter constant (Carter 1968). The above equation together with eqns. (j3J) 
and PJ provides a complete set of integrals of motion for the geodesic of a test particle in 
the Kerr geometry where Q, E and X z are three constants motion. It can be seen from the 
above equation that for a particle always confined to the equatorial plane (i.e., 9 = and 
9 = 7r/2), the Carter constant Q is zero. For a = 0, Q + \ 2 Z = A 2 is square of the total angular 
momentum of a particle's orbital motion. It is to be noted that the information about 9 and 
9 for an orbital trajectory is consolidated into a single parameter Q, which is a conserved 
quantity. Now using eqns. (j5J), ©, (j4]), (jSJ) and ([7]), we obtain the most general equations 
of motion governing the orbital trajectory of a test particle in the Kerr geometry (Bardeen 
1972; Wilkins 1972) given by 



= ±VR, (8) 
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= ~ aE + ~\ + + « 2 ) - (io) 

ar sin 6* A 

r/f r 2 -I- n" 1 

E-j- = -a 2 Esin 2 fl + aA 2 + ~T [E(r 2 + a 2 ) - A,a], (11) 
ar A 



where 



i? = [£(r 2 + a 2 ) - X z a} 2 - A [m 2 r 2 + (A, - aE) 2 + Q], (12) 

9 = Q-cos 2 #[a 2 (m 2 -£ 2 ) + (13) 

sin 

2.2. Effective potentials and prescription of orbital inclination 

R and are in general called the effective potentials describing the particle motion in 
a generalized Kerr spacetime. For a particle moving perpetually at the equatorial plane, 
the only effective potential that governs the particle motion is R with Q equals to zero. In 
order to understand the generic orbit, both potentials should be accounted for. To have an 
insight of the nature of orbits in threedimensional Kerr geometry, we recall that for bound 
geodesies (Q > 0) orbits either remain perpetually at the equatorial plane (Q = 0) or cross 
it repeatedly (Q > 0). More precisely, orbits with Q > (more general case), also called 
the off-equatorial orbits, are inclined to the equatorial orbital plane. Figure [T] shows that 
the orbital plane of an off-equatorial orbit is inclined at an angle % to the equatorial plane, 
which also means that the orbital plane makes an inclination % with the spin axis of the 
hole (symmetry axis). For a particular orbit, the inclination angle % is a constant of motion. 
However, in general, the plane of an off-equatorial orbit precesses with an angular frequency 
2a/r 3 called Lense-Thirring precession (Lense & Thirring 1918). Owing to this fact, the plane 
of the orbit neither remains fixed, nor crosses the equatorial plane with constant inclination 
i, but rather at a different angle. Looking meticulously, we find that the information of 
entire 9 motion (as for a particular off-equatorial orbit, 9 has all arbitrary values and it 
changes constantly) is glued up into one single inclination parameter i. Thus the motion of 
the bound orbit in the Kerr geometry is now governed by four parameters r, i, and t. 

As our intention is to derive a generalized PNP for particle orbits not confined to the 
equatorial plane, which itself is an approximate quantity, we restrict ourselves to an ideal 
case of circular orbits. In general, due to plane precession, the orbit never remains circular, 
but have eccentricity. The study of generic orbits (inclined and eccentric) for our particular 
case is deferred for future. The circular orbits satisfy R = dR/dr = 0. Hence to analyze the 
orbital motion for circular case, the radial effective potential R is necessary and sufficient. 
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Our next objective is to describe the inclination parameter i with known quantities. 
As stated earlier, the constant inclination angle % (without considering plane precession) is 
the angle between the plane of the orbit and the spin axis of the hole. As A 2 = Q + A 2 for 
a = stated in the preceding subsection, we identify in the Cartesian coordinate system 
Q ~ A 2 . + A 2 approximately. Hence Q can be treated as the component of square of the total 
angular momentum in perpendicular to the symmetry axis or to the spin axis (Aj_) of the 
black hole or precisely the component projected into the equatorial plane. However, even for 
a maximally spinning hole the error is found to be less than a few percent (Glampedakis et 
al. 2002). With this conjecture, the inclination parameter i can be defined as (Ryan 1995; 
Shapiro 1996; Hughes 2000) 



Taking all these into account, our radial effective potential R in eqn. fTTZj) with further 
simplification eventually reduces to 

R = r (r 3 + ra 2 + 2Ma 2 )E 2 - r(r - 2M)A 2 - 4aMr\ z E - A m 2 r 2 - A A 2 tan 2 i. (15) 

Thus the effective potential is a function of r and i (which is a conserved quantity for a 
particular orbit). If we consider the effect of Lense-Thirring precession on the orbital plane, 
then we need to amend a minor correction with the above definition of inclination angle. 
Now if you specify for circular orbits, then the effective potential further reduces as 



2aMA 2 + A 2 {r 2 



+ ^2sii(r 3 + ra 2 + 2Ma 2 )}A + m 2 r(r 3 + ra 2 + 2Ma 2 )A sec 2 i 



1/2 



(r 3 + ra 2 + 2Ma 2 ) sec % 



(16) 



which is essentially the solution for E of R = 0. 



2.3. Description of pseudo-Newtonian force 



The conditions for a circular orbit are 



R = 0, 



dR 

dr 



0. 



(17) 



Solving for E and A 2 from eqn. (TlTI) we obtain 




E 



A 




m 



m 
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E and X z are specific energy and specific angular momentum parallel to the symmetry 
axis respectively, and r is the radius of the circular inclined orbit. We do not represent 
the functions F\ and F 2 explicitly as they are cumbersome. It can be easily verified that 
A 2 sec i = A is the total angular momentum of the orbital plane with inclination angle i 
and its direction is perpendicular to the plane of the orbit. The relativistic specific angular 
momentum of the matter is defined by the relation t = —u^/ut = \ Z /E (Kozlowski et 
al. 1978). Following M02 and as standard practice, we can define the Keplerian angular 
momentum distribution A# = iseci. Therefore, corresponding centrifugal force in the Kerr 
geometry can be written as 

Xk 2 2MA 2 sec 2 i 



aV2VMr 3 / 2 {A + 2r(r - M)} + r A j [A + r 4 — a 2 (A + r 2 — 3Mr)] cos2i sec 2 i > 



1/2-, 2 



= F G = F Kr (19) 

where A = a 4 + r 4 + 2a 2 r(r — 2M). Thus from above, Fx r can be identified as the generalized 
gravitational force of the Kerr black hole at the Keplerian orbit. The subscript K represents a 
Keplerian orbit and subscript r represents force in the radial direction. The above expression 
reduces to that of M02 at an inclination angle % = 0, i.e. particle motion confined to the 
equatorial plane, and that of PW for a — i — 0. Thus we can say that above eqn. (fl9l) is 
the most general form of force corresponding to the PNP for a generalized accretion disk 
around a rotating black hole. One may be confused with the mathematical form of F^ r i n 
eqn. ffl9|) which appears complicated compared to that of PW and thus the PNP approach 
with Fjc r appears to be lacking the main motivation and advantage in comparison with 
full framework of GR. However, one should also remember that mathematical complexity 
of accretion disk dynamics in full Kerr geometry is beyond comparison with that in much 
simpler Schwarzschild geometry. Therefore, the apparent complex form of Fk t above is 
nothing compared to the complications exist in the set of actual equations in the Kerr 
geometry. Therefore, F^ r must be acceptable if it exhibits the essential features of the Kerr 
geometry what we examine in §3. 

The general form of corresponding PNP is 

V PK = [ F Kr dr. (20) 



Therefore, Vpx is a generalized axisymmetric potential mimicking the geometry of the general 
Kerr spacetime. The above potential can be attributed to the threedimensional free-fall 
acceleration of a test particle around a Kerr black hole. In Fig. [2] we show the evolution 
of the pseudo-Newtonian potential as a function of r and i for extremally rotating black 
holes. As the unstable region including unbound orbits with small inclination is broader for 
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a counter rotating black hole (described in detail in §3 with Tables 1 and 2), Vpk increases 
and tends to diverge at a larger radius than that of the corotating one. For the corotating 
case, Vpk tends to diverge at a very inner region in the vicinity of the black hole only. Figure 
[3] describes the corresponding contours of constant Vpk- It is found that the entire family of 
contour is shifted towards smaller r for the corotating case compared to the counter rotating 
one. Above qualitative features are same for other Kerr parameters. 

Now we write the Keplerian force components in the radial and vertical direction in 
cylindrical coordinate system (p, 0, z) from eqn. 0191) as 

2MA 2 B- 1 ' 2 

{A + p 4 B 2 - a 2 (A + p 2 B - UIpB 1 / 2 )} 1 ^] 2 

(21) 



aV^MpZ/W^A + 2pB 1 / 2 (pB 1 / 2 - M)} + A^p 2 - z 2 {A + p*B 2 - a 2 (A + p 2 B - ZMpB 1 ! 2 )} 



and 



2MA 2 B- l ' 2 z 

2 {A + p 4 B 2 -a 2 {A + p 2 B - ZMpB 1 / 2 )} 112 } 2 

" (22) 



^BV4{A + 2pB l / 2 { P B 1 l 2 - M)} + ApWy/fP- z 2 {A + p 4 B 2 - a 2 {A + p 2 B - 3MpB 1 / 2 )} 



respectively, where B = l+z 2 /p 2 , A = p 2 B-2M P B 1 / 2 +a 2 , A = a i +p i B 2 +2a 2 pB 1 / 2 (pB 1 / 2 - 
2M). Fk p and Fk z are the components of force corresponding to the pseudo-Newtonian 
vector potential components in cylindrical geometry. 

Once we have the pseudo- Newtonian potential (and force), we examine the correspond- 
ing density distribution d from the Poisson's equation 

1 d d 
-4nd = -— (pF Kp ) + —F Kz (23) 
pop oz 

shown in Fig. H]for extremally rotating black holes. It is found that the hypothetical density 
distribution d is Gaussian along z axis for both corotating and counterrotating black holes. 
Like Vpki d increases much faster at an outer radius for the counterrotating black hole than 
the corotating one. d for the corotating case increases significantly only at very close to the 
black hole event horizon. 



3. Comparison of the essential features for the Kerr geometry and 

pseudo-Kerr 

To check the validity of the potential derived in the preceding section, we follow the 
procedure considered in M02 and Mukhopadhyay & Misra (2003). We make a comparison 
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to show how far the potential Vpk reproduces the values of r s , and E s as that in Kerr 
geometry. We do not repeat the corresponding equations to calculate the parameters in Kerr 
geometry as they are established well in literature (see e.g. Shapiro & Teukolsky). However, 
for completeness, we show them for the pseudo-Kerr case. Stability of an orbit requires 
(PW) 

±(X K ) > 0. (24) 

Therefore, we obtain marginally stable orbit (r s ) by solving eqn. (|24p for r with "equals to" 
sign. 

The energy of a particle in the Kerr geometry according to Vpk is given by 

E= r -^ + V PK . (25) 
2 ar 

Therefore marginally bound orbit (rb) can be accounted for by solving the equation 

E\ n = (26) 

for rb- In addition, once we know r s , we can compute E\ rs = Es, the energy at last stable 
circular orbit. We enlist the comparison of r s , rb and E s between the Kerr and pseudo-Kerr 
geometry with percentage error for a few inclination parameters i in the Tables 1,2,3. We 
present the values of above parameters upto the orbital inclination i = 30°, as we believe 
that realistically the disk thickness does not exceed beyond the height corresponding to 
inclination parameter i ~ 30°. It is exciting to note that r s computed from pseudo-Kerr 
geometry exactly matches to that in Kerr geometry for all values of a and i. To best of our 
knowledge, none of the existing pseudo- Newtonian potentials proposed for threedimensional 
Kerr geometry reproduces r s exactly. Table 1 shows that values of obtained from Vpk 
match with that from exact Kerr geometry in very good agreement; with a maximum error 
~ 7% at high inclination angle 30°. With the increase of a upto certain value, error increases 
to a maximum (obtained in the range 0.998 > a > 0.9) and then decreases to at a = 1. It 
is to be noted that values of r s and r& are independent of i at a — 1 which coincide with the 
event horizon as the case for % = 0. The energy of an accreting particle at marginally stable 
orbit is also being reproduced by our PNP in very good agreement with that in exact Kerr 
geometry for % < 20°, as understood from Table 2. Although the error increases with the 
increase of i beyond 20°, except for a = 1 which is a special case, it is still within 20% for any 
realistic disk. For counter-rotating black holes the errors are less than that for co-rotating 
ones. However, it is to be noted that for very high inclination angle this trend reverses in 
E s as shown in Fig. [5] and Table 2. In both the cases of Kerr and pseudo-Kerr geometry, 
the orbits recede with the increase of inclination parameter i for prograde black holes, and 
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for retrograde black holes this dynamics is opposite. We compare the variation of specific 
energy for circular Keplerian orbit with the radial distance in Fig. [5] computed for exact 
Kerr geometry and our PNP for both prograde and retrograde holes at orbital inclination 
10° and 30°. 

The pseudo-Newtonian potential described by Semerak & Karas (1999), while tries to 
capture frame-dragging effects, can not reproduce other general relativistic properties e.g. 
radii of marginally stable and bound orbit, energy at marginally stable orbit. These are 
very important for the purpose of accretion disk producing jet, expected to be launched 
from inner edge of the disk. On the other hand, our potential is derived directly from the 
spacetime metric, more precisely from the effective potential of the threedimensional Kerr 
geometry given by eqns. (|15|) and ([TBI) . Either of the equations inherently carries information 
of frame-dragging effects with the term proportional to a\ z . Therefore, the frame-dragging 
effects are expected to be guaranteed by our potential. While the ad-hoc potential given by 
Chakrabarti & Mondal (2006) tries to mimic frame- dragging as well as inner disk properties 
both, it is not valid for entire parameter regime. It also does not reproduce the radii of 
marginally stable orbit exactly as ours does. While their potential is a function of r and 9 
both, they discuss most of the results at equatorial plane {6 = tt/2) and thus it is not clear 
how effective their potential is in a real system. Therefore, our potential is expected to be 
most selfcontained and useful to study generalized accretion disk including jet, particularly 
for the hydrodynamical/magnetohydrodynamical analysis. 

Thus, for a realistic accretion disk, our PNP, although approximate in nature, should 
reproduce the essential features of the spacetime in very good agreement. The potential is 
effectively useful to study the detailed hydrodynamics of the realistic disk with the inclusion 
of vertical outflow using Newtonian MHD equations, especially in the inner region, where 
general relativistic effects are very important. Thus the potential satisfies all the criteria to 
be claimed as good pseudo-Newtonian vector potential for a general Kerr metric. 

4. Discussion 

In this paper, we have prescribed a most generalized gravitational potential, pseudo- 
Newtonian vector potential, for the modeling of thick accretion disks around rotating black 
holes. As this potential is derived from the spacetime metric directly, it naturally carries 
several essential general relativistic properties. Our potential is valid for co-rotating as well 
as counter-rotating black holes with all realistic values of orbital inclination. The inherent 
nature of reproducing exact values of r s for the entire regime of Kerr parameter and for 
all inclination angle by our PNP makes it an efficient choice. The prescribed PNP also 
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reproduces in exact Kerr geometry with a maximum error ~ 7% for a high inclination 
~ 30°. However, it can be noted that the potential may not be a good approximation for 
the angular and epicyclic frequencies for less than 20 Schwarzschild radii. 

The full scale disk dynamics is very complicated and there is no virtue to assume that 
the accretion disk around black holes is fully Keplerian in nature. Moreover, the vertical 
outflow and wind in the form of relativistic jet, which are believed to be present almost in 
all quasars and microquasars, are predicted to be emanating from a very inner region of a 
puffed up disk or from a coronal transition region (Liu et al. 1999; Meyer et al. 2000). In 
addition, various inner fluid dynamical properties and also the powering of jet depend on 
spin parameter of the black hole. Thus evidently the realistic accretion physics (especially 
the accretion dynamics) around black holes can be studied with our prescribed axisymmetric 
generalized PNP, capturing the essential physics of GR. 

The general relativistic study of full scale dynamics of a realistic accretion disk for a 
highly turbulent viscous medium with the inclusion of jet and magnetic field is beyond the 
scope with the present physical and mathematical tools available to us. Most of the numer- 
ical studies of general relativistic flows have been performed with simplifying assumptions 
and remained confine to certain regions. Such simulations do not explicitly disclose the 
underlying microphysics of the structure and dynamics of disks and jets. In order to avoid 
complexities in full GR frameworks and with a notion to understand the underlying physics 
of the symbiotic connection of disk and jet, our PNP is highly recommended. Even if the 
results using our PNP deviate quantitatively by a few percent from that of the actual Kerr 
geometry, the accurate physical interpretation of the activities around black holes is possible. 
Our generalized PNP can be used to perform quasi-analytical studies and time-dependent 
simulations of accretion and wind and related physics like that done by Chan et al. (2005, 
2007), Lipunov & Gorbovskoy (2007), Shafee et al. (2007) for more realistic situations. 

Next, out of several plans, we would like to use our very approach to understand the 
physics of the symbiotic connection of disk and jet around black holes which, in the astrophys- 
ical context, is called by jet-disk symbiosis. In this context, one can extend the qualitative 
picture of the jet-disk symbiosis developed earlier (Falcke & Biermann 1995, Donea & Bier- 
mann 1996) to a full scale threedimensional MHD in the Newtonian framework with the 
inclusion of our PNP. The structure of the disk is highly modified in presence of the jet. 
Thus, we plan to describe the disk including all the above mentioned effects in a realistic 
situation, and to develop a tool to test with observation, both of the disk and the jet (from 
radio to TEV photons). 

The authors thank the referee for his/her very useful and constructive suggestions which 
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comments in preparing the final version of the manuscript. 
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Table 1: Values of r b 
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Note. — Percentage error=Er. Inclination angle i expressed in degrees, K and PK denote Kerr and 
pseudo-Newtonian result respectively. 
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Table 2: Values of E. 
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Note. — Percentage error=Er. Inclination angle i expressed in degrees, K and PK denote Kerr and 
pseudo-Newtonian result respectively. 
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Fig. 1. — Schematic diagram showing how an orbit of the accreting particle inclines to the 
equatorial plane. The orbit, inclined at an angle i, represents a generalized accretion. A 
represents the total angular momentum of the orbiting particle, a is the spin parameter of 
the black hole and z axis is the symmetry axis of the system. 
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Fig. 2. — Variation of generalized pseudo-Newtonian potential Vpk as a function of r (in 
units of GM/c 2 ) and i (in radian) for (a) a = 1, (b) a = —1. 
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Fig. 3. — Contours of constant Vpk for (a) a — 1, when solid lines correspond to Vpk = 

0.045, 0.06, 0.075, , 0.15, dotted lines to 0.2, 0.4, 0.6, , 1, and dashed lines to 1.5, 2, (b) 

a = —1, when solid lines correspond to V PK = 0.045,0.06,0.075, ,0.15, dotted lines to 

0.2, 0.4, 0.6, , 1, and dashed lines to 2, 4, 6, , 10. 
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Fig. 4. — Variation of density distribution generating Vpk as a function of r and z (in units 
of GM/c 2 ) for (a) a = 1, (b) a = -1. 
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Fig. 5. — Variation of specific energy for circular orbits as a function of radial coordinate 
at various Kerr parameters. The solid and dotted curves are for our pseudo-Newtonian 
potential with % = 10° and 30° respectively and the dashed and dot-dashed curves are for 
Kerr geometry with % = 10° and 30° respectively at each panel. 



